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Abstract 

We report a high-statistics measurement of the branching fraction for t~ ^ 'K~Tr^Vr and the 
invariant mass spectrum of the produced 7r~7r'' system using 72.2 fb~^ of data recorded with the 
Belle detector at the KEKB asymmetric-energy e^e~ collider. The branching fraction obtained 
is (25.24 lb 0.01 lb 0.39)%, where the first error is statistical and the second is systematic. The 
unfolded -k^tt^ mass spectrum is used to determine resonance parameters for the p{770), /o'(1450), 
and /?"(1700) mesons. We also use this spectrum to estimate the hadronic (27r) contribution to the 
anomalous magnetic moment of the muon (a^'^). Our result for a'^^ integrated over the mass range 
Vs = 2m^ - 1.8 GeV/c^ is aJJ" = (523.5 ib 1.5 (exp) ib 2.6 (Br) ib 2.5 (isospin)) x 10-^°, where the 
first error is due to the experimental uncertainties, the second is due to the uncertainties in the 
branching fractions and the third is due to the uncertainties in the isospin-violating corrections. 

PACS numbers: 13.40.Gp, 13.35.Dx, 14.60.Fg 
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I. INTRODUCTION 



Hadronic decays of the r lepton provide a clean environment for studying the dynamics of 
hadronic states with various quantum numbers. Among the decay channels of the r lepton, 
r~^7r~7r° has the largest branching fraction The decay is dominated by intermediate 
resonances and thus can be used to extract information on the properties of the p(770), 
p'(1450), and p"(1700) mesons and their mutual interference. 

From the conservation of vector current (CVC) theorem, the vr'vr*^ mass spectrum in 
this range can be related to the cross section for the process e"'"e" — tt+tt" and thus 
used to improve the theoretical error on the anomalous magnetic moment of the muon 

~ (S'm ~ 2)/2. Recent reviews of calculations of are given in Refs. [1, [s], 0]. It is 
known that the theoretical error on is dominated by the contribution from the leading- 
order hadronic vacuum polarization ajj'^*^'^°. This contribution cannot be derived within the 
framework of perturbative QCD and is usually evaluated using dispersion relations and the 
experimental cross section for e^e~ annihilation to hadrons [5|: 
relates the properties of the vr+vr" system produced in e^e~ - 



0,[i]. Alternatively, CVC 



TT ' TT 



to those of the vr 7r° 



system produced in r 



■TT TT 



Ur decay; thus, using CVC and correcting for isospin-violating 



effects, T data have also been used to obtain a more precise prediction for a^^""^'^^ 0, 0, i, Q ■ 
Recently, new precise data on e+e~ 7r"'"7r~ have become available from the CMD-2, 
KLOE, and SND experiments 0, [ll|, [H, [H, Q [H, 0, [IJ • ALEPH 
and OPAL [H, 

decay 



19 , CLEO 



21 



measured both the 27i spectral function and the branching fraction 
the latter was also determined by L3 2J] and DELPHI 



for the r — > TT 

Recent evaluations of the hadronic contribution to a 

-10 
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fj, usmg e~^e' 

while that using the r lepton data where applicable gives 
^ „^ (9.4 ± 10.5) X 10"^'' where the experimental value a^^P is dominated by the 
BNL E821 measurement [1^ (11 659 208.0 ± 6.3) x 10"^*^. These differences correspond to 
3.3 and 0.9 standard deviations, respectively. For the evaluation based on the e~^e~ data, 
a deviation of similar size that corresponds to a 3.4 cr discrepancy is claimed in Ref. (29| . 
To clarify these differences between the e"'"e~-based and r-based predictions, more data 
on e~^e~ — > ir'n'^ and r~ —>■ Tr'n^Ur decays are needed. In this paper we present a high- 



th 



(27.5 ± 8.4) X 10" 



data result in a°fP — a„ 



27 



a 



cxp 



a 



th 







statistics measurement of the n n mass spectrum produced in r 



■ 71 



decays using 



data collected with the Belle experiment at the KEKB asymmetric-energy e~^e~ collider 
operating at a center-of-mass (CM) energy of 10.6 GeV. The data sample is about 50 times 
larger than those of previous experiments. 



II. BASIC FORMULAS 



The differential decay rate for r ^ tt tt Vj- can be expressed as [30 



dr(r ^71 TT^ Ut) 

d^ 



-po ^^\^ud\ 'Jew 



— 

mi . 



2s 



mt 



with 



Glml 



(1) 



(2) 



" 1927r3 

Here s is the invariant mass squared of the 7r~7r° system, v_{s) is the weak spectral 
function characterizing the 7r~7r° system, Gp is the Fermi coupling constant, \Vud\ = 



4 



m-T 



1776.99t°:i 



0.97377 ± 0.00027 |31| is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element (32 



MeV/c^ |31| is the r lepton mass and 5"^^ accounts for short-distance 



electroweak radiative corrections for the vr vr^ system. The measured electron decay rate of 
the r lepton is related to Fg by 




r(r- ^ e-.- ^ r°5W = r° 1 + ^ ^ - , (3) 



where S^^^ is the electroweak radiative correction for the decay e~i7ei^T- 

The corresponding tt+tt" spectral function vo{s) can be obtained from the e~^e~' — > tt+tt" 
cross section 

a(e+e- ^ vr+vr-) = t;o(s), (4) 

where s is the e'^e^ CM energy squared and ao is the fine-structure constant at s = 0. Up 
to isospin-violating effects, CVC allows one to relate the spectral function from r decays to 



the isovector part of the e~^e spectral function [33 



v.{s) = vt\s). (5) 

Alternatively, the mass spectrum of the two-pion system can be expressed in terms of pion 
form factors, which are useful for comparing resonance shapes in the charged and neutral 
two-pion systems. The spectral function Vj{s) (j = — , 0) is related to the form factor F^{s) 
via 

vM = ^ms)\^ (6) 

where /3_(s) (/3o(s)) is the pion velocity in the vr'vr^ (vr+vr") rest system. The ve- 
locities Pj{s) are explicitly given by /3-(s) = A^''^(l, m^_/s,m^o/s) [3^ and Po{s) = 
X^/\l,ml./s,ml./s) = [l-4m2_/s], with X{x,y,z) = [x - + ^/^f][x - - ./^f]. 

The hadronic physics is contained within Vj{s) or, equivalently, in F^{s). One goal of this 
analysis is to provide a high-precision determination of the weak form factor using 
r~— > vr~vr° data, so that a comparison with from the e~^e~ data can be used to 

test CVC. From Eqs. ([1]), ([3]) and ([6]), one can obtain the basic formula that expresses the 
form factor F~ (s) in terms of the observables: 



where ;B,r7r is the branching fraction, (l/A'^7r)(dA'^^7r/ds) is the normalized invariant mass- 
squared distribution for the r^^vr~vr° Ur decay. Be is the branching fraction for e~Vri^e 
and S'ew = S^/S^^^. 

In this paper, we report new measurements for both the branching fraction iS^^r and 
the normalized mass spectrum (l/A^^^)(diV^^/ds). These results are used to provide a new 
evaluation of the hadronic contribution to the muon anomalous magnetic moment from the 
2vr channel. 
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III. DATA SAMPLE AND SELECTION CRITERIA 



The data sample used was collected with the Belle detector at the KEKB asymmetric- 
energy e^e~ collider j35|. It is based on an integrated luminosity of 72.2 fb^^ recorded at 
a CM energy of 10.58 GeV. The Belle detector is a large-solid-angle magnetic spectrometer 
consisting of several detector components. Charged track coordinates near the collision point 
are measured by a three-layer silicon-vertex detector (SVD) that surrounds a 2 cm radius 
beryllium beam pipe. Track trajectory coordinates are reconstructed in a 50-layer central 
drift chamber (CDC), and momentum measurements are made together with the SVD. An 
array of 1188 silica-aerogel Cherenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and specific ionization measurements {dE/dx) in the CDC 
provide a capability for the identification of charged particles. Photon detection and energy 
measurement of the photons and electrons are provided by an electromagnetic calorimeter 
(ECL) consisting of an array of 8736 CsI(Tl) crystals all pointing toward the interaction 
point. These detector components are located in a magnetic field of 1.5 T provided by 
a superconducting solenoid. An iron fiux-return located outside the coil is instrumented 
to identify muons and to detect mesons (KLM). A comprehensive description of the 



detector is given in Ref. [36 



To study backgrounds and determine selection criteria, we perform Monte Carlo (MC) 
simulation studies for various processes. Signal and background r"'"r~-pair events are sim- 
ulated using the KKMC generator 37|]. The r decays are modeled with the TAUOLA 
39| in which the values of the branching fractions are updated to more recent 
The cross section for e+e~ r+r~(7) is also updated to the recent measure- 




values [4 

ment reported in Ref. 



41 



by the PHOTOS program 



The radiative corrections to the r-hadronic decays are simulated 
4^ . The QQ generator 43] is used for BB and qq continuum 
processes, the BHLUMI 4^ program for radiative Bhabha events, the KKMC |37| program 
for radiative /i^/i^-pair events, and the AAFH 45|] program for two-photon processes. The 
BHLUMI and KKMC programs include higher-order radiative corrections and are among 
the most accurate programs available. The detector response is simulated by a GEANT3- 
based program [46|. In order to simulate beam-induced background realistically, detector 
hits taken from randomly triggered data are added to wire hits in the CDC and to energy 
deposits in the ECL. Uncertainties due to imperfections in the Monte Carlo generators and 
detector simulation are discussed in the later sections. 



A. T+r pair selection 

The event selection consists of two steps. Initially, a sample of generic e^e~ — r'^r~(7) 
events is selected with relatively loose criteria. From this sample r~ 7r~7r° Vj. decays 
are identified. The number of generic t^t~ events is used to determine the t" ^ tx^ti*^ Vr 
branching fraction. 

Generic r+r~ events are selected by requiring that the number of charged tracks in an 
event be two or four with zero net charge; that each track have a momentum transverse 
to the beam axis (pt) of greater than 0.1 GeV/c to avoid tracks reentering the CDC; and 
that each track extrapolate to the interaction point (IP) within ±1 cm transversely and 
within ±5 cm along the beam direction to suppress tracks that originate from beam-particle 
interactions with the residual gas in the vacuum chamber. To suppress background from 
Bhabha and fi~^fi~ events, the reconstructed CM energies and the sum of the momenta of 
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FIG. 1: Missing niciss (-^V^miss) versus the polar angle direction of tlie missing momentum (^^jgg) 
for (a) the data, (b) MC e^e~ — > t^t~ events, (c) MC Bhabha and ii^iJL~{'y) events and (d) 
two-photon processes. Events inside the octagonal region are selected as r'^r~-pair candidates. 



the first and the second highest momentum tracks are required to be less than 9.0 GeV/c. 
The maximum pt among the tracks is required to be greater than 0.5 GeV/c. Beam-related 
background is rejected by requiring that the position of the reconstructed event vertex be 
less than 0.5 cm from the IP in the transverse direction and less than 2.5 cm from the 
IP along the beam direction. The polar angle of the leading particle with respect to the 
beam axis {9*) in the CM frame is required to be in the fiducial region of the detector: 
35° <e* < 145°. 

To reduce the remaining background from Bhabha, {'-/), and two-photon events, 

a requirement is imposed in the plane of the missing mass Mmiss and the direction of 
missing momentum in CM 6'^;^^, where M^[ss is evaluated from the four- momenta of the 
measured tracks and photons: (Mmiss)^ = {Pm — Ptr — P'yY- In this expression pi^ is the four- 
momentum of the initial e~^e~ system, while ptr and p^ are the sum of the momenta of 
measured tracks and photons, respectively. A pion mass is assumed for the charged tracks 
if they are not identified as electrons or muons. Each photon (reconstructed from clusters 
in the calorimeter) must be separated from the nearest track projection by at least 20 cm 
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and have an energy greater than 0.05 GeV in the barrel region (—0.63 < cos^ < 0.85), and 
greater than 0.1 GeV in the endcap region (—0.90 < cos^ < —0.62 and 0.85 < cos^ < 0.95). 
Photons near the edge of the detector fiducial volume are rejected. Scatter plots of Mmiss 
versus ^miss for data, the r-signal MC, the Bhabha and the two-photon MC are shown in 
Figs. [I]-(a), (b), (c) and (d), respectively. The band of events in data at M^iss ~ is due 
to backgrounds from Bhabha and fi~^fi~{'y) processes. Small vertical bands at 0^^^^ = 45° 
and = 150° are Bhabha events where the energy of one of the final state electron/positron 
is poorly measured because it has scattered in the material at the boundary of the barrel 
and endcap calorimeters. The events in the high-Mmiss region (> 7GeV/c^) are from the 
two-photon processes. 

Events within the octagonal region are selected as r"''r~ candidates to avoid the tail from 
background processes. 

^10 

EM 

o 
> 



(A 
(0 

c 





50 100 150 

6*miss (degree) 

FIG. 2: Missing iTiciss (il^fmiss) versus the polar a-iiglc direction of the missing momentum (^j^jgg) for 
data. The solid (dashed) lines show three vertical (horizontal) slices that are used to present the pro- 
jections in Fig. [3l The coordinates of the vertical (horizontal) lines are ^miss ~ ^0°' ^5°' 125°, 160° 
{M^iss = 0.5 GeV/c2, 2.0 GeV/c^, 5.0 GeV/c^, 7.0 GeV/c^). 

To display the r-pair and background contribution quantitatively, we divide the scatter 
plots of Mjniss vs ^miss iiito three vertical and three horizontal slices as shown in Fig|2l 
Projections for the six slices are shown in Fig. [3l where each process shows a characteristic 
shape: the t~^t~ candidates dominate in the central region in Mmiss and ^^iss- Both Bhabha 
and fi~^fi~ show a prominent peak at Mmiss ^ 0, but the width for the Bhabha is slightly 
wider than that of /U^/x". We use the events in the region |Mmiss| < 0.5 GeV/c^ to determine 
the normalization for the Bhabha and /i"'"/^". As two-photon processes dominate in the high- 
Mmiss region, the normalization for the two-photon processes is determined using the events 
at |Mjniss| > 8.0 GeV/c^. The arrows with solid (dotted) lines indicate the narrowest (widest) 
areas used to select r+r~-pairs by the octagonal selection. Although overah features of the 
data are modeled reasonably well by MC, some discrepancies are seen, for example, in 
the regions Mmiss ~ and > 6 GeV/c^ in Fig|3]-(c), which are taken into account as the 
systematic error on the background estimation. 
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FIG. 3: Projections to the missing mass (Mmiss) and the missing direction (^miss)- corre- 
spond to the vertical shces from left to right in Figl2j (d)-(f) correspond to the horizontal slices 
from top to bottom. The solid circles represent the data, and the histogram represents MC simu- 
lation (signal+ background). The open histogram shows the contribution from r"'"r~-pair process, 
the vertical (horizontal) striped area shows that from two-photon leptonic (hadronic) processes; 
the wide (narrow) hatched area shows that from Bhabha {^^fi~) process; and the shaded area 
shows that from the qq continuum processes. The arrows with solid (dotted) lines indicates the 
widest (narrowest) region corresponding to the octagonal boundary shown in Fig[2j 
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Acoplanarity 



FIG. 4: Characteristic distributions for surviving t'^t^ candidates: (a) Mmiss, (b) ^miss, (c) particle 
multiplicity Xpart = ('^tr + '1-7)1 X {n^r + n^)2-, (d) acoplanarity angle ^. The points indicate the 
data, the open histogram shows the r-pair MC and the hatched histogram shows the background 
from e'^e~ — > qq and other sources. All selection criteria are applied for (a) and (b). All criteria 
except for the quantity in question are applied for (c) and (d). The arrows in (c) and (d) indicate 
the boundary used to select a r-pair sample. 



Candidate events are divided into two hemispheres in the CM frame by the plane per- 
pendicular to the highest momentum particle, and the remaining background from e^e~ 
annihilation is suppressed by selecting events with low multiplicity as characterized by the 
quantity Xpart = {n^^ + ^^7)1 x ('^tr + ^^7)25 where ritrj- and the numbers of tracks 

and photons in hemisphere j. We require Xpart < 25. Finally, in order to eliminate Bhabha 
events in which one or both electrons produce a shower in material near the interaction 
region, the acoplanarity angle ^ between the first and second highest momentum tracks is 
required to be ^ > 1°, where ^ = | |0i — 02| — 7r| is defined as the two-track acoUinearity in 
azimuth. The Xpart and ^ distributions after applying all selection criteria except for the 
quantity in question are shown in Figs. |4]-(c) and (d), respectively. The selection boundary 
is shown by the arrows. 
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-8 -6 -4 -2 2 4 6 



(m,„ - m o)/a 

FIG. 5: Normalized 77 invariant mass (S^^) spectrum for data (points) and the r~— > h~TT^ Vr signal 
MC (open histogram), for the sample described in the text. The data plotted here correspond to 6% 
of the full data sample used in this analysis. The arrows indicate the signal region — 6 < S-y-y < 5 
and the sideband regions 7 < < 9. The sideband regions are used to subtract fake-vr'^ 

background. The shaded histogram shows the non-r background determined from MC simulation. 



After applying all selection criteria, 22.83 x 10^ r'^r^-pairs survive. The Mmiss and 6'miss 
distributions for the surviving events, shown in FigJ3]-(a) and (b), respectively, demonstrate 
a low level of the background and an overall good agreement between the data and the MC 
model. 

For surviving events, the dominant background is from the e^e~ qq {q = u, d, s, c) 
continuum and amounts to (5.30 ± 0.53)% of the total number of events. The systematic 
error for the qq background is determined from the uncertainty of the normalization of the 
events in the region 25 < Xpart < 30, where qq processes dominate. The background from 
e'^e~ — > T(4S') — > BB is small (0.1 %). Backgrounds from Bhabha, two-photon 
leptonic and hadronic events are to be 0.92 ± 0.09 %, 0.28 ± 0.01%, 0.62 ± 0.03% and 
0.60 ± 0.09%, respectively. Here the systematic errors for each background is determined 
from the uncertainty of the normalization of the events in the background enhanced region 
mentioned before. 

B. —>■ ir^ TT^ Uj- selection 

Within the r"''r~-pair sample, ^ ti^tt^ Ur decays are reconstructed by requiring that 
there be both one charged track and one ir^ in a single hemisphere. The vr" candidate is 
selected based on the normalized invariant mass 5*^^ = (m^^ — m^o)/cr^-y, where a^^ is the 
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FIG. 6: Invariant-mass-squared (M^^o) distribution for t~— >7r~7r'^ Vj- after imposing tight tag-side 
requirements. The sohd circles with error bars represent the data, and the histogram represents 
the MC simulation (signal -|- background) . The open histogram shows the contribution from r~— >■ 
7r~7r''p'T-; the narrow cross-hatched area shows that from t~ ^ K~ ir^ 1/^; the wide cross-hatched 
area shows that from r~ — > h~ {mr^)^^; and the striped area shows that from the qq continuum and 
other non-r processes. 



mass resolution of the 77 system. The value of ranges from 0.005 GeV/c^ to 0.008 
GeV/c^, depending on the vr*^ momentum and polar angle. Pairs of photons with jS'^^l < 9 
are considered as tt" candidates. To keep beam-related background at a negligible level, we 
require that the CM momentum of the vr^ be greater than 0.25 GeV/c and the photon CM 
energy be greater than 0.08 GeV. 

The distribution of S^^ for the selected tt^tt^ sample, with one charged track and one 
TT^ candidate in a single hemisphere, is shown in Fig. [51 The lower-side tail of the S^.y 
distribution is primarily due to rear and transverse leakage of electromagnetic showers out 
of the CsI(Tl) crystals and the conversion of photons in the material located in front of 
the crystals. Good agreement between data (points) and MC (open histogram) indicates 
that these effects are properly modeled by the MC simulation. We define the interval 
—6 < 5*^^ < 5 as the vrO si gnal region. Spurious 7r° background is small and estimated from 
the sideband regions 7 < \S^^\ < 9. To reduce feed-down background from multi-Tr*^ decays 
such as T~ 7r~{n7r^)i>r { n > 2), signal candidates are rejected if there are additional 7's 
in the same hemisphere with energy greater than 0.2 GeV. 

The 7r~7r'' invariant-mass-squared (M^^o) spectrum is obtained assuming the pion mass 
for the charged track; it is shown in Fig. [6] along with the MC prediction. To improve 
the vr" energy resolution, a ir^ mass constraint is imposed. The spurious 7r° background 
level depends on the M^^o region, varying from 4% to 7%. (This is subtracted using 
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sidebands.) The final sample contains 5.43 x 10^ r~ h~Ti^ candidates after the 7r° 
background subtraction, where hr denotes vr" or K~ . This sample is 50 times larger than 
those of previous studies. 

The spectrum is dominated by the p(770) peak and a shoulder due to the p'(1450). A 
small but clear structure from the p"(1700) is visible at M^^o ~ 2.7 (GeV/c^)^. 

There are two sources of background: feed-down from other r decay modes and the 
gg-continuum. Feed-down background arises mainly from multi-vr^ modes such as r~ — >■ 
h-{mr'^)i^r (6.02 ±0.08%), r A'L/i-7r°i^^(0.48 ± 0.04%) and r ^ uj7i^Ur{uj 7r%) (0.10 ± 
0.01%). Here h~ denotes either vr" or K~ . After all modes are included, the total feed- 
down background level is (7.02 ± 0.08)%. The error given here includes a MC statistical 
uncertainty as well as the uncertainty on relevant branching fractions. The contribution of 
these feed-down backgrounds dominates at low values of M^^o (Fig- E]). 

The gg-continuum background level is (2.22 ±0.05)% in total, and is concentrated mostly 
in the high M^^o region above 2.0 (GeV/c^)^. Since the reduction of this high-mass back- 
ground is essential in the measurement of the mass spectrum, we impose the stringent 
requirement that the tag side contain only one charged track and no photons. This require- 
ment improves the signal-to-noise ratio in the high-mass region M^^o > 2.0 (GeV/c^)^ by a 
factor of 3, although the total size of the r"— s>7r~7r° Ur sample is reduced by a factor of 2.5. 
The normalization of the continuum MC is validated using data in the mass region above 
the r lepton mass: M^^o > i^r- Background from the other non-r processes, such as BB, 
Bhabha and /i^yU~7 in the final sample is negligible (< 0.1%). 



IV. MEASUREMENT OF THE BRANCHING FRACTION 
A. Basic Method 

The branching fraction for T~—>-h~7T^ {Bh^o) is determined by dividing the signal yield 
^hnO by the total number of selected r leptons 2A^^t- taking into account various efficiencies 
and background corrections: 

N^^O (1 _ fefecd-down _ ^non-r) / X ^ 

BhTrO = X ^ ; ^ X X -YTT ■ (8) 

In this formula, is the background fraction in the r+r" sample, e,-^ is the efficiency of 
the r"'"r~-pair selection, e^^o is the efficiency for r~ h^n^u decays to pass the r+r~-pair 
selection, and e]^o is the efficiency for t~ ^h'TT^v decays satisfying the r"'~r~-pair selection 
to pass the h~Ti^ selection. The product e^^o ■ e\^o is the overall detection efficiency for the 
h^-K^v final state. The parameter ft^ced-down ^j^g fraction of h^Ti^v candidates coming from 
other r decay modes, and h^°'^"^ is the fraction coming from non-r processes. In this formula, 
several common uncertainties such as that on the luminosity, on the cross section for t^t~- 
pair production, on the trigger efficiency, and on the t^t~ selection efficiency cancel in the 
ratio. In the measurement of the branching fraction, the stringent tag-side condition is not 
imposed to avoid any possible bias that it might introduce. The values for all factors are 
listed in Table [T] along with the MC statistical error. 
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TABLE I: Values of parameters used for the branching fraction measurement along with MC 
statistical errors. 



Parameter 


Value 


Err 


32.59 ± 0.05 % 


Jb = 


36.24 ± 0.07 % 
1.112 ±0.003 




41.01 ±0.13 % 


brr 


7.80 ± 0.03 % 


j,fccd-down 


7.02 ± 0.08 % 


Lnon-T 


2.22 ± 0.05 % 



B. Systematic uncertainties 

The sources of systematic uncertainties on Bh-KO are listed in Table [III The uncertainty 
on the tracking efficiency is estimated using D*^ — > D^n^ — >■ K~'k^'k~^ decays to be 1% per 
track. A large part of this uncertainty cancels in the ratio of Eq. ([H]); the resulting relative 
uncertainty from this source is IS.B/B = 0.47 %. 

TABLE II: Systematic uncertainties for the T~^h' branching fraction. 



Source of uncertainty 


Afi,,.o (%) 


i^B/B) (%) 


Tracking efficiency 


0.12 


0.47 


vr*^ efficiency 


0.32 


1.27 


Background for t~^t~ 


0.15 


0.59 


Feed-down background for r~ 


^h-^'K^Vr 0.04 


0.16 


Non-r background for r~— > h" 


7r° 0.05 


0.20 


7 veto 


0.05 


0.20 


Trigger 


0.08 


0.32 


MC statistics 


0.02 


0.08 


Total 


0.39 


1.52 



The systematic error on the 7r° detection efficiency has two components: one is the uncer- 
tainty coming from the 7r° selection criteria and the other is that from the absolute efficiency 
calibration. For the uncertainty coming from the vr" selection, we check the uncertainty by 
changing the definition of the signal and background region, by taking into account the 
uncertainty in the resolution function and by changing the 7r° threshold momentum. For 
example, the relative branching fraction changes by only ^B/B = 0.1% if the signal region 
is changed from the nominal one to —7 < S^^ < 7. Also the uncertainty is AB/B = ±0.2% 
for the changes of the 7r° threshold by ±0.05 GeV from the nominal value. 
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In order to make an absolute efficiency calibration independently of the signal process, 
we use the 77 — >■ 77 and rj — > tt^tt^tt'^ signals, whose branching fractions are known rather 



precisely. Combining the PDG world average 3l| for the 77 — >■ 77 and 77 —» vr'^vr^vr^ branching 



fractions and the recent measurement from the CLEO collaboration [47| , we obtain the ratio 
of the branching fractions of 



B{rj TiHH^) 



0.829 ±0.007, 



which has 0.84% relative accuracy. 

By comparing the signal ratio Ri = N{ri -k^-k^-k^ / N{ri 77)) for the data [i = 1) and 
the MC (z = 2), the correction factor for the detection efficiency of one ir^, rjcor, is determined 



to be ?7cor = y-RdataZ-RMC = 0.950 ± 0.012, where the error includes the uncertainties in the 
rj signal measurement (1.2%) and the errors on the rj decay branching fractions (0.4%). 

This correction factor is also con&med by a study of electron/positron tracks from photon 
conversions (i.e. 7— s>e"^e~) in the SVD region. It is found that the E/P distribution for 
those tracks is simulated correctly above 1.0 GeV, but requires some tuning below 1.0 GeV. 
This imperfection of the MC primarily leads to a difference in the 7r° signal shape and an 
efficiency difference between data and MC. 

The non-r background is dominated by qq continuum processes; this is estimated by 
using the events above the r mass: M^^o > fn^. The statistics of the data and MC sample 
determine the error. 

The uncertainty on the feed-down background /^^fced-down (.Qj^gg f^om the MC statistics 
and the uncertainty on the branching fractions for r~ h~{n7T^)i'r, t~ —>■ K~7r^Ur and 
^ um^Ur (cj ^ 7r°7). 

The veto of additional 7's is required in the event selection to reduce background from 
multi-TT*^ decay channels. However, this veto can reject signal itself if photons are radiated 
in the initial or final state and those photons are detected within the detector fiducial 
volume. In addition, photon candidates can also appear due to electromagnetic shower 
fragments and/or misreconstructed electrons. Therefore a precise simulation of the photon 
radiation as well as the shower simulation are important. The uncertainty from these sources 
is estimated by changing the veto threshold by ±0.1 GeV around the nominal value of 
0.2 GeV; the resulting relative change in i3^^o is only ±0.20%. Signal events are flagged by 
several trigger conditions that require two or more CDC tracks with associated TOF hits, 
ECL clusters, or a significant sum of energy in the ECL. This redundancy allows one to 
monitor the efficiency of each trigger requirement. The uncertainty arising from the trigger 
is estimated by assuming that there is a ±3% uncertainty on the track and energy trigger 
efficiencies, which is the maximum variation measured during experimental running. The 
resulting relative uncertainty is small (0.32%) since the t~^t~ trigger efficiency is high (97%). 



C. Results 

Inserting all values into Eq. ([8]) we obtain 

Bh^o = (25.67 ± 0.01 ± 0.39)%, (9) 

where the ffist error is statistical and the second is systematic. This result is in good 
agreement with previous measurements, as shown in Table IIIII Our statistical error is 
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significantly lower than those of the other measurements; our systematic error is similar to 
those of CLEO, L3 and OPAL, and larger than those of ALEPH and DELPHI. 



TABLE IIL Branching fractions for r ^h' measured by different experiments. 



Experiment 


BhA%) 


Reference 


CLEO 


25.87 ±0.12 ±0.42 




L3 


25.05 ±0.35 ±0.50 


m 


OPAL 


25.89 ±0.17 ±0.29 


[23] 


ALEPH 


25.924 ± 0.097 ± 0.085 


[19] 


DELPHI 


25.740 ±0.201 ±0.138 


[25] 


This work 


25.67 ±0.01 ±0.39 





We combine the PDG world average for the t~ ^ K'h^Vt branching fraction 3l| with a 
recent BaBar measurement 48|] to obtain the result Bx-nO = (0.428 ± 0.015)%. Subtracting 
this from our r~— >-/i~7r° Ur result gives a r~— >7r~7r'' branching fraction of 

B^^o = (25.24 ± 0.01 ± 0.39)%, (10) 

which is consistent with the previous measurements from CLEO 21 1 and ALEPH lol . 



V. MEASUREMENT OF THE MASS SPECTRUM 

In order to obtain the true 7r~7r'' mass spectrum, one must apply corrections for: (1) back- 
ground, (2) smearing due to finite resolution and radiative effects, and (3) mass-dependent 
acceptance. 



A. Background Correction 

As noted earlier, there are three sources of the background that enter the r~^7r~7r° Ur 
sample: (1) fake vr" background, (2) feed-down background from other r decay channels, and 
(3) the background from the qq continuum. The total magnitude of these background con- 
tributions is about 7% in the p(770)) peak region, but the fraction of the background varies 
strongly with M^^o] there is approximately a 4-order-of- magnitude difference between the 
signal level in the p(770) peak region and that in the high M^^o region above 2.5 (GeV/c^)^. 
Thus a reliable estimation of the background is important for the measurement of the mass 
spectrum. 

The sidebands of the M^^ distribution are used to estimate the fake vr^ contribution. This 
background dominates at values of M^^o less than about 0.25 (GeV/c^)^. 

In the 7r~7r'' system, the feed-down background dominates at similarly low values of M^^o 
while the gg-continuum background dominates at high values of M^^o (see Fig. E]). These 
backgrounds are subtracted bin-by-bin. 
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FIG. 7: The acceptance determined from MC simulation as a function of the vr vr mass squared. 



N 3.5 




FIG. 8: Correlation between the generated and observed invariant masses squared of the vr vr' 
system in the t~— >vr~vr'^ Ut decay. 



B. Acceptance Correction 

The acceptance determined from MC simulation as a function of the generated vr^vr^ 
mass squared is shown in Fig. [71 The acceptance varies smoothly and its average value is 
7%. This acceptance includes a factor for the tag-side branching fractions B{t~ — > £~i/^z/^) 
(£~ = ^~,e~), B{t" h~UT-) {h~ = '7i~,K~), which does not affect the shape of the mass 
spectrum. The acceptance decreases at low values of M^^^ due to the overlap of 7 clusters 
with the vr" track in the calorimeter. 
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The detector effects include M^^o-dependent acceptance and bin-by-bin migration caused 
by tlie finite mass resolution. The radiative decay r~ — > ir'n'^'yi'r also causes some bin 
migration. We correct for these effects using an unfolding procedure that makes use of the 
MC to characterize the acceptance and the bin migration. These effects can be characterized 
by the acceptance matrix A defined by 

b = Ax, 

where x is the vector containing the generated 7r~7r'^(7) mass-squared spectrum and b is the 
reconstructed one. It is possible to apply the inverse of A to the spectrum observed in the 
data to obtain an unfolded spectrum. However, this procedure is not robust with respect to 
the statistical fluctuations entering the determination of A, and can yield unphysically large 
point-by-point fluctuations. To cure this problem, we use an unfolding program employed in 
the ALEPH experiment [49|. In this program, the unfolding is based on the Singular- Value- 
Decomposition (SVD) method, in which the acceptance matrix is inverted by constraining 
the number of singular values to only those elements that are statistically significant. 

The acceptance matrix is determined iteratively using a signal MC based on the 
KKMC/TAUOLA program. In the second iteration, the p"(1700) resonance is included 
in the MC based on our measurement. Final state radiation in r hadronic decays is simu- 
lated by the PHOTOS program. In order to take into account the effects of 7 radiation in 
the decay -k^tx^^Vt-, the invariant mass squared of the 7r~7r°7 system is taken as the 

generated quantity. 

The output of the program is the unfolded distribution and its covariance matrix. The 
correlation between the generated quantity and the measured one is shown in Fig. [HI The 
figure shows a clear correlation between the measured and generated values. The resolution 
in Ml^o is 0.005 (GeV/c^)^ in the low-mass region and 0.030 (GeV/c^)^ in the high-mass 
region; thus by choosing the bin size to be AM^ = 0.050 (GeV/c^)^, the off-diagonal com- 
ponents of the acceptance matrix are small. 



C. Systematic Uncertainties 

The sources of systematic errors associated with the unfolded mass spectrum 
(l/A^)(dn/ds) are subdivided into several classes according to their origin, which are the 
unfolding procedure (UNF), the background subtraction (BKG), the acceptance correction 
(ACC), and the energy scale (ENG). These contributions are summarized in Table HVl for 
each Min region and are described below. 

The systematic error due to the unfolding procedure is determined from MC by comparing 
the true and the unfolded results (UNFl). Another estimate of the uncertainty of the 
unfolding is made by changing the value of the unfolding parameter that determines the 
optimum number of the singular values of the acceptance matrix (UNF2). 

The uncertainty of the background subtraction is estimated for each source. BKGl is 
from continuum processes. Its uncertainty is estimated using the control sample in the mass 
region higher than the r mass. The statistics of the data and MC sample determine its 
error. BKG2 is the feed-down background. Its uncertainty is estimated by varying the 



branching fraction values 3l| used in the MC by ±lo". BKG3 is the non-7r° background. 
The uncertainty of the non-vr'' background is estimated by changing the vr^ sideband region. 
This uncertainty dominates in the threshold region but is negligible elsewhere. 
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TABLE IV: Relative systematic errors (in %) of the unfolded spectrum for each M^^o region for the 
different sources of uncertainty: unfolding procedure (UNFl, UNF2), the background subtraction 
(BKGl, BKG2, BKG3), the acceptance correction (ACC), and the photon energy scale (PES). See 
the text for a more detailed description. 



M a region 


l^irst bm 
(0.08) 


rni 1 1 J 

i hresnold 
region 
(0.2-0.3) 


P 

region 
(0.55-0.60) 


(1 0-1 2) 


/ 

P 

region 
(1 9-2 0) 


P 

region 
(2.5-2.7) 


UNFl 


2.50 


0.79 


0.31 


0.85 


1.50 


1.50 


UNF2 


2.60 


0.53 


0.09 


0.27 


0.58 


9.19 


BKGl 


1.13 


0.09 


0.01 


0.04 


0.52 


5.76 


BKG2 


4.90 


0.65 


0.10 


0.10 




0.50 


BKG3 


25.21 


4.80 










ACC 


5.36 


1.44 


0.03 


0.15 


0.15 


0.40 


PES 


1.24 


1.08 


0.59 


0.99 


0.05 


0.50 


Total 


26.5 


5.3 


0.7 


1.5 


1.8 


11.4 



The acceptance uncertainty is dominated by the uncertainty of the tt^ efficiency. This is 
estimated by changing the measured values of the photon efficiency by 1 standard deviation. 
In addition, the effect of requiring that photons be isolated from charged tracks is checked 
by changing the isolation criteria from 20 cm (default) to 30 cm. 

The uncertainty of the photon energy scale (PES) is estimated from the vr" peak position 
to be ±0.2%. This uncertainty is important for the peak position of the resonances. The 
uncertainty in the charged track momentum scale is negligible compared to that of the 
photon energy scale. 

Individual components of the uncertainty are added in quadrature to obtain total sys- 
tematic errors of 5.3% in the threshold region, 0.7% near the p(770) peak and 1.8% in the 
vicinity of the p"(1450) (see Table IV]). 

D. Results 

The unfolded s = M^^^o^^f ) spectrum dNj^T^/ds is shown in Fig. [9l The error bars in the 
figure include both statistical and systematic errors added in quadrature and in most cases 
they are smaller than the size of the data points shown by closed circles. The results are also 
presented in terms of the normalized unfolded spectrum (l/7V^^)(dA''^^/ds) in Table IVl and 
in terms of the pion form factor in Table IVI[ In these tables, the statistical and systematic 
errors are given separately. The statistical errors in the figure and the table are the square 
roots of the diagonal components of the covariance matrix. 

In Fig. [9l the p peak and a shoulder due to the p'(1450) are clearly visible. The dip at 
s ~ 2.5 (GeV/c^)^ is caused by destructive interference between the p'(1450) and p"(1700) 
resonances. 

To determine the parameters of the p, p' and p resonances, a fit using Breit-Wigner 
(BW) functions is performed. The pion form factor is parametrized with Breit-Wigner 
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FIG. 9: The fully corrected M^^o distribution for r~ ^ tt^tt" z^,-. The solid curve is the result 
of a fit to the Gounaris-Sakurai model with the p(770), p'(1450), and p"(1700) resonances. All 
resonance parameters (mass, width, phase and the normalization factor 1^7^(0)^) are allowed to 
float. 



functions corresponding to the p, p'(1450), and p"(1700) resonances: 

1 



FJs] 



1 + /5 + 7 



pi 1 



(11) 



where the parameters (3 and 7 (denoting the relative magnitude of the two resonances) are in 
general complex. We use the Gounaris-Sakurai (GS) model [sOl for the Breit-Wigner shape: 



BW, 



GS 



M? + d- M.TAs] 



(Mf-s) + f{s)-t^Us)' 



with an energy-dependent width 



r.(.) = r.f^^ 



(12) 



(13) 



Here, k{s) = ^y/sP-{s) is the pion momentum in the n rest frame. The functions /(s) 
and h{s) are defined as 



M} 



k\s) {h{s) - h{M^)) + (M^ - s)k\M^) 



ds 



s=M? 



(14) 
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TABLE V: The unfolded normalized spectrum (l/A'7r7r)(dA'^jr7r/ds) as a function of the invariant 
mass squared s = M^^q. The square of the diagonal element of the error matrix is taken for the 
statistical errors. Note that the scale is different for the left-and right-sides. 



Bin 


M2 


1 dN 
N ds 


Stat. 


Syst. 


Bin 


M2 


1 dN 
N ds 


Stat. 


Syst. 






X 10-^ 


X 10-=* 


X 10-^ 


No. 




X 10-* 


X 10-* 


X 10-* 


No. 


(GcV/c2)2 






GcV-2c4 


No. 


(GeV/c2)2 


GcV-2c4 


GeV-^c* 


GeV-2c* 


1 


0.088 


8.1 


3.1 


2.1 


32 


1.625 


341.47 


5.91 


5.83 


2 


0.125 


45.8 


2.4 


9.6 


33 


1.675 


290.94 


5.56 


5.39 


3 


0.175 


108.4 


2.4 


15.5 


34 


1.725 


250.39 


4.99 


4.32 


4 


0.225 


185.1 


2.6 


14.3 


35 


1.775 


210.22 


4.86 


3.40 


5 


0.275 


278.0 


2.8 


8.3 


36 


1.825 


170.63 


4.36 


2.75 


6 


0.325 


396.0 


3.3 


4.1 


37 


1.875 


139.72 


3.89 


2.26 


7 


0.375 


628.6 


4.2 


5.0 


38 


1.925 


109.26 


3.82 


1.85 


8 


0.425 


1024.2 


5.5 


14.7 


39 


1.975 


81.85 


3.20 


1.39 


9 


0.475 


1710.8 


7.2 


21.5 


40 


2.025 


63.08 


2.84 


1.25 


10 


0.525 


2643.0 


9.0 


22.5 


41 


2.075 


45.02 


2.72 


0.89 


11 


0.575 


3268.0 


9.7 


21.4 


42 


2.125 


29.89 


2.22 


0.65 


12 


0.625 


2755.5 


9.0 


19.4 


43 


2.175 


20.06 


1.94 


0.61 


13 


0.675 


1907.2 


7.4 


17.0 


44 


2.225 


13.08 


1.68 


0.42 


14 


0.725 


1214.0 


5.8 


8.4 


45 


2.275 


7.93 


1.25 


0.42 


15 


0.775 


826.4 


4.6 


10.7 


46 


2.325 


4.85 


1.82 


0.51 


16 


0.825 


592.8 


3.6 


8.6 


47 


2.375 


3.04 


2.84 


1.39 


17 


0.875 


435.6 


2.9 


3.5 


48 


2.425 


2.32 


3.36 


1.17 


18 


0.925 


327.5 


2.4 


4.6 


49 


2.475 


2.09 


2.39 


0.76 


19 


0.975 


253.7 


2.1 


3.7 


50 


2.525 


2.07 


1.14 


0.37 


20 


1.025 


206.3 


1.8 


2.7 


51 


2.575 


2.24 


0.64 


0.19 


21 


1.075 


172.5 


1.5 


3.6 


52 


2.625 


2.88 


0.78 


0.17 


22 


1.125 


144.0 


1.3 


1.4 


53 


2.675 


3.22 


0.94 


0.33 


23 


1.175 


124.0 


1.2 


3.7 


54 


2.725 


3.52 


0.95 


0.51 


24 


1.225 


108.7 


1.1 


3.3 


55 


2.775 


3.49 


0.89 


0.57 


25 


1.275 


94.0 


1.0 


1.0 


56 


2.825 


3.29 


1.60 


0.62 


26 


1.325 


82.0 


0.9 


2.4 


57 


2.875 


3.09 


1.86 


0.80 


27 


1.375 


72.6 


0.9 


1.3 


58 


2.925 


2.65 


1.81 


0.65 


28 


1.425 


64.5 


0.8 


1.8 


59 


2.975 


1.79 


0.79 


0.43 


29 


1.475 


56.1 


0.8 


1.8 


60 


3.025 


0.88 


0.39 


0.22 


30 


1.525 


48.0 


0.7 


1.4 


61 


3.075 


0.43 


0.20 


0.12 


31 


1.575 


41.1 


0.6 


0.7 


62 


3.125 


0.09 


0.04 


0.03 
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TABLE VI: The form factor squared as a function of the invariant mass squared s. The 

results are obtained by inserting the measured value (l/A^7r7r)(dA'^7r7r/ds) into Eq. d?]). Note that 
the short-distance radiative correction is already applied, where the value of Sew is taken to be 
1.0235 lb 0.0003 (see the discussion in the Appendix). 



Bin 


M2 


IF P 


Stat. 


Syst. 


Bin 


M2 


IF P 

1 TT 1 


Stat. 


Syst. 


No. 










No. 


(GcV/c^f 








1 


0.088 


1.434 


0.549 


0.377 


32 


1.625 


0.711 


0.012 


0.012 


2 


0.125 


1.707 


0.091 


0.358 


33 


1.675 


0.636 


0.012 


0.012 


3 


0.175 


2.362 


0.053 


0.337 


34 


1.725 


0.576 


0.011 


0.010 


4 


0.225 


3.211 


0.045 


0.248 


35 


1.775 


0.511 


0.012 


0.008 


5 


0.275 


4.260 


0.042 


0.127 


36 


1.825 


0.439 


0.011 


0.007 


6 


0.325 


5.622 


0.046 


0.058 


37 


1.875 


0.382 


0.011 


0.006 


7 


0.375 


8.492 


0.057 


0.067 


38 


1.925 


0.318 


0.011 


0.005 


8 


0.425 


13.392 


0.072 


0.193 


39 


1.975 


0.255 


0.010 


0.004 


9 


0.475 


21.894 


0.093 


0.275 


40 


2.025 


0.211 


0.009 


0.004 


10 


0.525 


33.384 


0.113 


0.284 


41 


2.075 


0.162 


0.010 


0.003 


11 


0.575 


40.996 


0.122 


0.269 


42 


2.125 


0.117 


0.009 


0.003 


12 


0.625 


34.503 


0.112 


0.243 


43 


2.175 


0.085 


0.008 


0.003 


13 


0.675 


23.936 


0.093 


0.214 


44 


2.225 


0.060 


0.008 


0.002 


14 


0.725 


15.324 


0.074 


0.106 


45 


2.275 


0.040 


0.006 


0.002 


15 


0.775 


10.525 


0.058 


0.137 


46 


2.325 


0.027 


0.010 


0.003 


16 


0.825 


7.637 


0.047 


0.111 


47 


2.375 


0.019 


0.018 


0.009 


17 


0.875 


5.693 


0.038 


0.046 


48 


2.425 


0.017 


0.024 


0.008 


18 


0.925 


4.350 


0.032 


0.061 


49 


2.475 


0.017 


0.019 


0.006 


19 


0.975 


3.435 


0.028 


0.050 


50 


2.525 


0.019 


0.011 


0.003 


20 


1.025 


2.851 


0.024 


0.037 


51 


2.575 


0.024 


0.007 


0.002 


21 


1.075 


2.439 


0.022 


0.052 


52 


2.625 


0.036 


0.010 


0.002 


22 


1.125 


2.087 


0.019 


0.020 


53 


2.675 


0.050 


0.014 


0.005 


23 


1.175 


1.847 


0.018 


0.056 


54 


2.725 


0.066 


0.018 


0.010 


24 


1.225 


1.667 


0.017 


0.050 


55 


2.775 


0.083 


0.021 


0.013 


25 


1.275 


1.486 


0.016 


0.015 


56 


2.825 


0.102 


0.050 


0.019 


26 


1.325 


1.339 


0.015 


0.039 


57 


2.875 


0.132 


0.079 


0.034 


27 


1.375 


1.229 


0.015 


0.022 


58 


2.925 


0.165 


0.113 


0.040 


28 


1.425 


1.132 


0.014 


0.031 


59 


2.975 


0.178 


0.079 


0.043 


29 


1.475 


1.025 


0.014 


0.032 


60 


3.025 


0.165 


0.073 


0.041 


30 


1.525 


0.913 


0.013 


0.027 


61 


3.075 


0.203 


0.098 


0.056 


31 


1.575 


0.818 


0.013 


0.014 


62 


3.125 


0.287 


0.138 


0.081 
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with dh/ds\j^2 = h{Mf) \{8k\Mf)y^ - (2Af2)-ij + (27rAf2)-i and 



^ 3 ml / M + 2fc(M/) \ _ m^M, 

7rF(M2) 1^ 2m, j 2'nk{Mf) iik^Mf)' ^ ' 

Note that the function d is chosen so that the BW^^ function is unity at s = [501] . 

Since the unfolded mass spectrum has bin-by-bin correlations, the off-diagonal compo- 
nents of the covariance matrix X are included in the evaluation: 

= E - /(^^; «)) {y, - f{s,- a)) , (17) 

where Ui is the measured value at the i-th bin, f{s]a) is the value of the function for 
parameters a, and {X~^)ij is the inverse of the covariance matrix. 

There are 10 parameters in this formula: the masses (Mj) and widths (Fj) for the p, p', 
and p" resonances, their relative amplitudes I7I, and phases and (p^. In addition, as 
an overall normalization factor, we introduce |F7r(0)p as an additional parameter. In the 
BW form, this value should be unity. However, in order to take into account a possible 
deviation from the form of the fitting function, two kinds of fits, in which this parameter is 
either fixed or floated, are carried out. The other 10 parameters are floated in the fit. 

The results of the fit are shown as the solid line in Fig. [9] for the M^^o distribution as well 
as in Fig. [inland Fig.fTlTa) and (b), where the results are compared directly to the weak form 
factor squared \F~{s)\^, derived bin-by-bin from Eq. ([7|). In Eq. ([?[), we use the world average 
value (including our measurement) for the branching fraction i^jrvr = (25.24 ± 0.10)% and 
for the CKM matrix element V^d = 0.97377 ± 0.00027 31i- For the short- distance radiative 



correction Sew = S^w/^ew^ ^^^^ value 1.0235 ± 0.0003, to be consistent with the 
isospin breaking correction discussed in Ref. jl, [13] (see the appendix for more details). 

The fitted results are summarized in Table [ylTI for the cases when |F,(0)p is fixed to 
unity (the second column) and is allowed to float (the third column). In the table, the first 
error is statistical and the second is systematic. The value of the P^r degree of freedom 
(NDF) is 80/52 for the fixed and 65/51 for the floated cases. It is found that |-F,(0)p is close 
to unity (|F,(0)p = 1.02 ± 0.01 ± 0.04) even when it is allowed to float. It should be noted 
that the data can be fitted using BW resonances only, without any additional background 
terms. The fit quality for the fixed case is slightly worse than that for the floated one. The 
curves shown in Figs. [9l [TO] and [TlT a) andfTlTb). correspond to the case when the parameter 
1^7^(0) p is floated, but the differences between the floated and the fixed cases are small. 

The significance of the p"(1700) signal is given in the last row of Table WUl The signifi- 
cance is determined from the change in the when the signal and its associated degree of 
freedom are removed from the fit. If the p"(1700) signal is excluded from the fit, the fo'^ 
the fit increases by 55 (60) units, in the case that |F(0)p is fixed to unity (allowed to float). 
This increase in the for the fit, with the joint estimation of four removed parameters 
(mass, width, I7I, 0^), corresponds to a 6.5cr (7.0cr) significance for the p"(1700) signal (sij . 
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FIG. 10: Pion form factor for t~^tt~tt^ Ur- The solid circles are the Belle result while the squares 
and stars show the result of ALEPH [l9| and CLEO 2C ] , respectively. The error bars for the Belle 
data include both statistical and systematic errors added in quadrature. The solid curve is the 
result of a fit to the Gounaris-Sakurai model with the /3(770), /9'(1450), and p"(1700) resonances, 
where all parameters are floated. 
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FIG. 11: (a) Pion form factor for r~^7r~7r'^ z/^ in the p(770) mass region and (b) in the threshold 
region. The convention of the plots is the same as in Fig. [TOl The error bars for the Belle data 
include both statistical and systematic errors added in quadrature. The solid curve is the result of 
a fit to the Gounaris-Sakurai model, where all parameters are floated. See the text for details. 
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TABLE VII: Results of fitting the „ distribution for r — > vr vr'' z^T- to the Gounaris-Sakurai 
model with the p(770), /9'(1450), and /?"(1700) resonances. The results are shown for two cases, 
fixed |-F,r(0)P = 1 (the second column) and all parameters are allowed to float (the third column). 
The first error is statistical and the second one is systematic. The systematic errors include the 
uncertainty in the backgrounds, unfolding as well as the uncertainty of the photon energy scale. 
The last row gives the significance of the p"(1700) signal. 



Parameter 


Fit result 


Fit result 




ffixed iFfOlPl 


1 all free 1 

I ail iL y 


Mp, MeV/c 


774.6 ± 0.2 lb 0.5 


774.9 ± 0.3 ± 0.5 


ip, MeV 


148.1 ± 0.4 ± 1.7 


148.6 ± 0.5 ± 1.7 


M„i MeV/c^ 


1446 ± 7 ± 28 


1428 ± 15 ± 26 


Tp,, MeV 


434 ± 16 ± 60 


413 ±12 ±57 


4)13, degree 


0.15±0.05l[j:J^ 
202 ± Af^l 


0.13±0.0llg:J^ 

197 ± 91:^^ 


Mpn, MeV/c2 


1728 ± 17 ± 89 


1694 ± 41 ± 89 


Tp,,, MeV 
l7l 

(/)^, degree 


164 ± 2ll^^ 
0.037 ± O.OOetaooi 
24 ± 9i^28 


135 ± Sdtle 
0.028 ± 0.020+[j:|]^^ 
-3 ± 13l^i 


|F(0)p 


[1.0] 


1.02 ±0.01 ±0.04 


xVndf 


80/52 


65/51 


p"(1700) signif., a 


6.5 


7.0 



The fit parameters are correlated, with the correlation matrix: 
/ 



1.00 


















0.58 


1.00 
















0.45 


0.39 


1.00 














-0.30 


-0.13 


0.24 


1.00 












-0.14 


-0.28 


-0.15 


0.36 


1.00 










-0.31 


-0.24 


-0.06 


0.47 


0.13 


1.00 








-0.42 


-0.08 


0.25 


0.34 


-0.61 


0.38 


1.00 






0.31 


0.22 


0.08 


-0.10 


-0.41 


0.37 


0.29 


1.00 




0.32 


0.27 


0.13 


-0.13 


-0.44 


0.48 


0.28 


0.64 


1.00 


0.54. 


0.29 


-0.01 


-0.07 


0.05 


0.52 


-0.27 


0.54 


0.67 


-0.19 


-0.06 


0.13 


0.61 


0.23 


0.49 


0.57 


0.53 


0.17 



where the parameters are |F„ 



(18) 



M, 



pi 



Tp, Mpi, Fp/, |/?|, 0^, Mp//, Fp//, I7I and 



The systematic uncertainties for the parameters are summarized in Table IVIIIl where the 
sources of the systematic errors from the unfolding procedure (UNF), the qq background 
subtraction (BKGl), the feed-down background subtraction (BKG2), the acceptance cor- 
rection (ACC) and the photon energy scale (PES) are shown separately. The uncertainty 
in the p mass (0.7 MeV/c^) is mainly due to the uncertainty in the photon energy scale. 
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TABLE VIII: Systematic errors for resonance parameters from the unfolding procedure (UNF), 
the background subtraction (BKG), the acceptance correction (ACC) and the photon energy scale 
(PES). 



Mp Tp Mp, Tp, (t)p Mpn Tpu \^\ cj)^ 

(MeV/c2) (MeV) (MeV/c^) (MeV) (deg.) (MeV/c^) (MeV) (deg.) 



Fit Bias 


0.3 


1.6 


25 


49 


0.028 


4 


75 


10 


0.006 


13 


UNF 

BKGl 

BKG2 


0.3 
0.3 


0.3 


4 
11 
1 


24 
25 


0.020 

+0.143 
-0.031 


4 

+41 
- 5 

2 


11 
13 
2 


14 

+86 
-10 

2 


0.002 

+0.053 
-0.020 

0.001 


12 

+117 
- 22 


ACC 




0.1 


1 


4 








7 




1 


PES 


0.3 


0.6 


2 


1 




2 


45 


15 




1 


Total 


0.5 


1.7 


28 


60 


+0.147 
-0.047 


+41 
- 8 


89 


+89 
-26 


+0.053 
-0.021 


+118 
- 28 



The uncertainty in the qq background (BKGl) dominates for the p'(1450) and p"(1700) 
resonance parameters. 

The values of the p(770) mass and width are consistent with the results of the previous 
measurements in r decay. For the p'(1450), a slightly higher mass value than in the previous 
measurement is obtained. It is found that this value is sensitive to the value of the resonance 
parameters for the p"(1700) since they both interfere. This is the first time that all the 
parameters for the p(770), p'(1450) and p"(1700) are determined in a single fit. Production 
of the p"(1700) in r~ decays has been unambiguously demonstrated and its parameters 
determined. 



E. Comparison of Belle and previous r data 

Comparisons of the pion form factor squared |F^(s) p measured in Belle to those measured 
by CLEO HqI and ALEPH [l^ experiments are given in Fig. [TlTa) for the p(770) region and 



in Fig.fTlTb) for the low-mass region M^^o < 0.4(GeV/c^)^. Figure [T2] shows a more detailed 
comparison, the difference in |F^(s)p of Belle and CLEO, ALEPH data for the fit of the 
Belle data divided by the fit value, for the mass-squared range 0.11-1.20 (GeV/c^)^. It can 
be seen that over the entire mass range shown, the |F^(s)p values from Belle are consistent 
within errors with those of CLEO. Agreement is worse when Belle data are compared to 
ALEPH data. Below 0.5 (GeV/c^)^ our points are mostly higher than those of ALEPH, 
while above 0.7 (GeV/c^)^ they are systematically and significantly lower. 



VI. IMPLICATIONS FOR THE MUON ANOMALOUS MAGNETIC MOMENT 
A. Basic Formula 

As described in the introduction, the hadronic vacuum polarization term aha,d,LO p^g^yg qji 
important role in the standard model prediction for the muon anomalous magnetic moment 
ttp] the error on this contribution is the most significant source of the uncertainty in a p. In 
this section, we discuss the implication of our measurement for the determination of a^^'^'^'^. 
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-0.2 



-0.3 




FIG. 12: Comparison of the pion form factor squared IFttI-s)]^ measured with the Belle detector to 



results from the CLEO [20(] and ALEPH 19| experiments in the p{770) and p'(1450) mass region. 
The difference from the fit of the Belle r~^7r~7r^ Uj- data divided by the fit value is plotted. The 



solid circles show the Belle result, the open stars show the CLEO result [20|] and the hatched band 



19[. Both Belle and ALEPH results include the systematic errors. The 



shows the ALEPH result 
inner error bars in the Belle data indicate the statistical errors. 




The value of a^^'^'^'~' is related to the e^e annihilation cross section via the dispersion 
integral 

(19) 

where a^^^ is the total cross section for e~^e~ — >hadrons at the center-of-mass energy ^/s. The 
superscript in a^^^ denotes the "bare" hadronic cross section, which is defined as the mea- 
sured cross section, corrected for QED radiative corrections such as initial state radiation, 
electron- vertex correction and the vacuum polarization in the photon propagator. However, 
the final state radiation (FSR) photons coming from the process cr°(e^e~ — > V'y — > 7r"'"7r~7) 
should be included in the a^^^ term (See the detailed discussion on page 12of Ref. 0). The 
kernel function K{s) is given by 

K(s) = 
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with X 



4m2/ s . K{s) is a smooth function increasing 



from 0.63 at the threshold s = 4m^ to unity at s = oo. Because of the 1/s dependence of 
'^hadi^) additional 1/s factor in the integral in Eq. (IT^ . low-mass hadronic final 

states dominate the contribution to a^^'^'^^: in fact about 70% of a^atixo jg ^^^g ^j^g 
pion final state with 4m^ < s < 0.8 (GeV/c^)^. Consequently, the 2n mass spectrum from 
r data is useful for obtaining predictions for a^ad.LO ygjj^g CVC. 



TABLE IX: Summary of the aJJ" contribution from r data from Belle, ALEPH [JJ] and CLEO [2C 
experiments. The errors are only shown for the Belle data. In this table, the isospin breaking 
correction is not applied except for the short-distance radiative correction given by the term Sew- 
The errors in the Belle data are only statistical. The errors for other experiments are similar or 
slightly worse than the Belle value in the same mass range. 



ranee 








(10-10) 




Integrated 




10) 


((GeV/c2)2) 




Belle 


ALEPH 


CLEO(*) 




Belle 


ALEPH 


CLEO(*) 


0.075-0.200 


39.55 




0.97 


38.20 


43.81 


39.55 




0.97 


38.20 


43.81 


0.200-0.350 


70.62 




0.46 


66.84 


75.71 


110.3 




1.07 


105.0 


119.5 


0.350-0.500 


123.25 




0.28 


119.10 


106.10 


233.5 




1.11 


224.1 


225.6 


0.500-0.650 


196.78 




0.23 


194.00 


197.80 


430.3 




1.13 


418.2 


423.3 


0.650-0.800 


62.35 




0.10 


62.35 


73.03 


492.6 




1.14 


480.5 


496.3 


0.800-0.950 


15.64 




0.04 


16.40 


14.31 


508.3 




1.14 


496.9 


510.6 


0.950-1.100 


5.74 




0.02 


6.50 


6.01 


514.0 




1.14 


503.4 


516.6 


1.100-1.250 


2.86 




0.01 


3.27 


2.32 


516.9 




1.14 


506.7 


519.0 


1.250-1.400 


1.65 




0.01 


1.89 


2.65 


518.5 




1.14 


508.6 


521.6 


1.400-1.550 


1.03 




0.01 


1.13 


0.77 


519.5 




1.14 


509.7 


522.4 


1.550-1.700 


0.60 




0.00 


0.66 


0.49 


520.2 




1.14 


510.3 


522.9 


1.700-1.850 


0.36 




0.00 


0.37 


0.33 


520.5 




1.14 


510.7 


523.2 


1.850-2.000 


0.19 




0.00 


0.19 


0.21 


520.7 




1.14 


510.9 


523.4 


2.000-2.600 


0.13 




0.00 


0.21 


0.15 


521,8 




1.14 


511.1 


523.6 


2.600-3.200 


0.13 




0.00 


0.21 


0.04 


522.0 




1.14 


511.1 


523.6 



(*) For the CLEO 



data, the boundary of the mass range is slightly different from the one shown in 
the first column, since ^/s bins are used. 



B. Results 

Details of our determination of aj^^, the 27r contribution to a^^'^'^'~' , are given in the 
appendix, where the basic formulas, the corrections applied for the isospin-violating effects, 
and discussions on the error estimation are presented. 

Our result on a'^^ over the mass range y/s = 2771^^ — 1.8 GeV/c^ is 

a7[2m^, 1.8 GeV/c^] = (523.5 ± 1.5 (exp.) ± 2.6 (Br.) ± 2.5 (isospin)) x 10"!°, 
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where the first error is due to the experimental uncertainties, i.e. the statistical error 
(1.1 X 10~^°) and experimental systematic error (1.0 x lO"^'') added in quadrature. The 
second error comes from the uncertainties in the branching fractions. The third one is 
the error on the isospin-violating corrections. These sources of error are discussed in the 
appendix. 

This result can be compared to those from previous ALEPH, CLEO, and OPAL r data. 
The combined result given in Ref. is 

al''[2m^, 1.80 GeV/c^] = (520.1 ± 2.4 (exp.) ± 2.7(Br.) ± 2.5 (isospin)) x 10"^° (r) . 

In terms of the experimental error (i.e. the first uncertainty), our result improves the 
previous combined result by 40%. A detailed comparison from our results and those of 
ALEPH ^ and CLEO is given in Table HX] for the a)^^ contribution. As seen in Table HXl 
the contribution from the mass-squared region 0.8 (GeV/c^)^ < s < 1.25 (GeV/c^)^, where 
one observes a deviation between the Belle and ALEPH data in Fig. 9, is only 4.6% of the 
total 2tt contribution. In this region, the contribution from the ALEPH measurements is 
higher than that from Belle but it is compensated by the opposite tendency in the region 
s < 0.50 (GeV/c^)^. Moreover, since the contribution from CLEO is between Belle and 
ALEPH, the difference between our result and the combined previous-r result becomes 
smaller. Consequently, our result agrees well with that of the combined result given from 
the previous r data within < la of the experimental error. 

On the other hand, the value of aJJ^ in the same ^/s region evaluated from the e~^e~ cross 



section measurements is 27 



a7[2m,, 1.80 GeV/c'] = (504.6 ± 3.1 (exp.) ± 0.9 (rad.)) x lO"''' (e+e" : CMD2, SND) , 

where the first error includes both statistical and experimental systematic errors added in 
quadrature. The second error is due to radiative corrections. 

Our r result is noticeably higher than the e^e~ result. This confirms the longstanding 
difference between the spectral functions of the 2-7r systems produced in r-decay and e+e~ 
annihilation 5^ • 

In summary, we have studied the decay r~— s>7r~7r° using a high-statistics data sample 
taken with the Belle detector at the KEKB e"'~e~ collider. The branching fraction is mea- 
sured with 1.5% accuracy. From Table UTTl we can calculate the accuracy of the previous 
experiments: CLEO 1.7%, L3 2.4%, OPAL 1.3%, ALEPH 0.5%, DELPHI 0.9%. These com- 
parison shows that the accuracy of the Belle result is better than CLEO and L3, similar to 
OPAL and worse than ALEPH and DELPHI. The result is in good agreement with previous 
measurements. In the unfolded 7r~7r° mass spectrum, in addition to the p(770) and p'(1450) 
mesons, the production of the p"(1700) in r~ decays has been unambiguously demonstrated 
and its parameters determined. The unfolded spectrum is used to evaluate the 2tt contribu- 
tion to the muon anomalous magnetic moment aj^^ in the region y/s = 2m.,^ — 1.80 GeV/c^. 
Our results agree well with the previous r based results but are higher than those from e^e~ 
annihilation. 
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Appendix 



C. Determination of a^,'^ from Belle data 

In this appendix, we list the details of our determination of a^^ over the range -y/s = 
2m^ — 1.8 GeV/c^. We discuss the basic formula, the integration procedure, the corrections 
for the isospin-violating effects, and the evaluation of the experimental errors. 



1. Basic formula 

In the isospin symmetry limit, CVC relates the quantities in r^— > tt^tt'^ z/^ decay to the 



cross section e+e tt+tt (ctJJ^) through the relation [30 



^TTTrlCVC 



M{s)n 



1 fB^ 

X 



ds 



1 diV^. 



(20) 



where N'{s) is given by 



mud? 



mi 



2s 



(21) 



A more precise link between hadronic spectral functions from r decays and the e^e~ 
hadronic cross section requires a calculation of radiative corrections as well as the inclusion 
of the isospin breaking effects (both of kinematic and dynamic origin). These effects have 
recently been discussed by Cirigliano et al. 33] and by Flores-Tlalpa et al. [s^j- According 
to them, the formula in Eq. fl2Ul) is modified to 



I CVC 



N{s) 



X 



1 diV^, 
iV^TT ds 



'Rm{s) 
S- 



EW 



(22) 



with 



Rm{s) 



F,(s) 



(23) 



where Sew is the dominant short- distance electroweak correction and Rib{s) takes care 
of the isospin-violating corrections. Rib{s) includes the long-distance QED correction 
Gem{s), the phase space correction factor (3q{s)/ (3^{s), and the ratio of the pion form factors 
\Fo{s)/F4s)\\ 



2. Integration procedure 

Using the measured distribution (l/A^7r7r)(dA''^7r/ds), the moment a'^^ can be obtained by 
inserting the bare cross section Eq. ( !22l) to Eq. ( |T9l) and integrating over s. The integration 
in Eq. (fT9l) is carried out numerically by taking the sum of the integrand evaluated at 
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the center of each bin. The statistical error on a^J^ is calculated including the off-diagonal 
elements of the covariance error matrix X, 



E 



dai 



(24) 



There are several external parameters in these equations; the values used for them are 
listed in Table 1X1 For m^, Vud, and Be, PDG [3l| values are used. For the ir'n^ branching 
fraction, our measurement is consistent with the world average given in Ref. 3l|. Including 
our result and the recent ALEPH i3^^o measurement [l^, the new world average is 



(25.42 ±0.10)%. 



(25) 



We use this new world average for the evaluation of aj^^. 

The errors on a'^^ arising from external parameters are summarized in Table |Xl the total 
systematic error from these sources is ±2.7 x 10"^'^ (dominated by ABt^j^o). 



TABLE X: Values of the external parameters and systematic errors on aJJ^ arising from these 
sources. 



Source 


Value 


Relative error 
(%) 


(10-10) 


Reference 


Vud 


0.97377 ± 0.00027 


0.027 


± 0.26 


m 


Be 


(17.84 ± 0.05)% 


0.28 


± 1.45 


m 


Btttt" 


(25.42 ± 0.10)% 


0.41 


± 2.13 




Total external 






± 2.6 





3. Experimental systematic uncertainty 

The systematic errors on a^'^ arising from internal sources (specific to this measurement) 
are listed in Table IXII and discussed below. There are two sources of background in the 
7r-7r° sample: (i) feed-down from r" — ^ h~{n7r^)h'T-, r" — ^ K'tt^Vt, t~ — >■ ^^'^^{uj — >■ n^'^) 
and (ii) non-r background. In the first case, MC statistics and the uncertainty on the 
branching fraction are used to estimate the error. In the second case, the uncertainty on the 
background as estimated from the control samples is assigned as an error. As mentioned 
earlier, the fake-vr^ background is subtracted using sideband events and the uncertainty is 
determined by varying the signal and sideband regions. 

It is found that the shape of the mass spectrum is insensitive to uncertainties in the 7r° 
efficiency, as it is only at a few % level. The uncertainty of the integration procedure comes 
from the binning effects. Adding all individual errors in quadrature we obtain a total error 
on a^^ arising from internal sources of ±1.0 x 10"i°. 

To check the stability of a^'^, we perform the following tests: 

1. The sample is divided into subsamples based on the tag-side topology, i.e., one electron, 
one-prong, or three-prong. The values of a^J^ obtained from these subsamples are 
consistent within the statistical errors. 
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TABLE XI: Systematic errors on d^^ arising from internal sources (specific to this measurement). 



Source 




Background: 




non-T (e+e" qq) 


± 0.11 


feed-down h{n-K^)v 


± 0.09 


feed-down K^-n^v 


± 0.15 


Energy scale 


± 0.10 


7r''/7 selection 


± 0.24 


7 veto 


± 0.93 


Efficiency: 




vr77 


± 0.35 


charged track 


< 0.10 


Integration procedure 


< 0.10 


Total internal 


± 1.04 



2. The sample is divided into subsamples based on the running period. Again, the values 
of a"^ obtained are consistent within the statistical errors. 

3. The sample might be sensitive to the requirement on the overlap region between the 
charged track and 7 clusters. To estimate this sensitivity, we select events with a 
tighter isolation requirement on 7's and on the track extrapolation: 30 cm instead of 
20 cm. 

The resulting variation in a"^ is small and is included as an additional systematic 
error. 



TTTT 



Jf.. Iso spin- violating corrections for a 

Three identifiable sources of isospin breaking corrections are the mass difference of the 
charged and neutral pions, p — u interference effects and the radiative corrections, which are 
included in the factor -Rib (s) /Sew in Eq. (120|1 . The size of the isospin- violating corrections 
from these sources and the possible uncertainties from other sources are summarized in 
Table IXlIl 

(i) The dominant contribution from electroweak radiative corrections comes from the 
short-distance corrections. In the leading logarithmic order, the short- distance radia- 



tive corrections to the decays r (du)^^- are enhanced by the factor [5J], [55 



S(m^, mz) = l + ^"^"""V + 2Q) In ^ = 1.01878, (26) 
Air mi 



where mz is the Z boson mass, and a{mr) = 1/133.50(2) is the QED coupling at the 

'2 _ In 

^3 3' 



r lepton mass. Q is the average quark-doublet charge. Therefore, Q = |(| — -1) — ^ 
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TABLE XII: Sources of the isospin violation between the e'^e and r spectral function in the 2tt 



channel and the corrections to a^^ . 



The correction is based on the procedure given in Refs. 30. l53l| . 



Source 


of 


Correction to Uncertainty on 


References 


isospin 


violation 


^nn (10-10) aJJ" (10-1°) 





Short-distance rad. cor. (Sew) 


- 12.0 


± 0.2 


Long-distance rad. cor. {Gem{s)) 


- 1.0 




m^- / m^o in phase space) 


- 7.0 




p — LO interference 


+ 3.5 


± 0.6 


m^- 7^ m^o (/?^ in the decay width) 


+ 4.2 




Electromagnetic decay modes 


- 1.4 


± 1.4 


vripO 7^ nip- 




± 2.0 


Total 


- 13.7 


± 2.5 



[5^, [55], [56j, [5j 
[30], [53] 



for the semileptonic decays, r Urud. Since Q = — | for leptons, there are no leading 
logarithmic corrections for leptonic decays. 

We can go further and sum up all short-distance logarithms of the a" In" niz via the 
renormalization group. This procedure replaces Eq. fl2^ by 



S{mr,mz) 



a[mz) 
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3 0!(mT) 
25 TV 




3 

23 IT 


a^(m^)J 









L01907. 



(27) 



where the last two terms are the short- distance QCD corrections |56| . 

Taking into account the subleading correction for the leptonic decay, the short- distance 
Electroweak correction Sew is given by 

1 



Sew = S{mr,mz) ^ 



L0235± 0.0003, 



(28) 



1 CIS, sub 

where b^^^ 



(l + "^^^-^ — TT^^^ = 0.9957. The difference between the resummed 
value (1.01907) and the lowest-order estimate (1.01878) is taken as the error on Sew 
and the resummed Sew value is used throughout this paper. The shift of a^'^ from 
this correction is Aa^'^ = (12.0 ± 0.2) x 10-^'^. Note that this correction is already 
applied in Eq. ([7j) and the pion form factor given in Table IVI[ 



(ii) The long-distance QED radiative correction Gem{s) was computed in the framework of 
chiral perturbation theory by Cirigliano et al. in Ref. [l^] • Recentl y, it was reevaluated 
based on a meson dominance model by Flores-Tlalpa et al. in Ref. |53| . It is found that 
the predictions of both models coincide if the contribution of the a;(782) intermediate 
state r ujti~'Vt-{uj vr'^7), is excluded in the latter model calculation. Since in our 
data this intermediate u;(782) contribution is already subtracted at the analysis level, 
we use the s-dependent correction factor provided by Flores-Tlalpa, which does not 
include the u) contribution. This correction produces a shift of Aa^^^ = — 1.0 x 10" 
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(iii) The vr and vr*^ mass difference in the ratio of the phase space (3q / /3i results in a shift 



of Aar = -7.0 X 10"^° , 



(iv) The 71 and vr" mass difference also affects the (3^ factor in the energy-dependent decay 
width (Eq. [T3D, which provides a positive shift Aa^'^ = +4.2 x 10~^° . 

(v) The effect of the p — u interference is estimated using the interference amplitude 
parameterized in the following form 30 



F'As),-. = ■ (29) 

3m^ m;f, - s - imj:^ 

For the numerical evaluation, we take dp^ = (—3.3 ± 0.7) x 10"'^ GeV^, = 0.783 



GeV and = 0.00844 GeV. The net effect is Aa^ = (+3.5 ± 0.6) x 10~^° , 



(vi) The largest source of the uncertainty for the isospin-violating effects is from the p~ 
and mass difference. The mass difference is consistent with zero within about 1 
MeV, which gives an uncertainty on aj^^ of ±2.0 x 10"^''. 

(vii) Finally, electromagnetic decay of the p meson is the source of the isospin violation. The 
decay irir'-f deserves particular attention. The shift and its uncertainty are estimated to 
be Aaji^'^ = (-1.4 + 1.4) x 10"^° from the width difference r(p° n+n-j) - T{p+ 
7r+7r°7) = (0.45 + 0.45) MeV [sj. 



Summing all these corrections, the overall isospin-violating correction and its uncertainty 
are estimated to be (—13.7 ± 2.5) x 10~^°. These corrections relate the r spectral function 
to the "pure" e+e~ tt+tt^ with all QED corrections switched off. Since the a^^^ term in 
Eq. fll9l) must include processes with the final state radiation (FSR) photons coming from 
the process e~^e~ — Vj — > 7r"''7r~7, we must reintroduce the FSR contributions j^. After 
including those contributions (+4.2 x 10"^'') [5|, the total correction becomes (—9.5 + 2.5) x 
10-10. 

With the external parameters in Table |X] and corrections discussed above, we obtain a'^^ 
over the range ^/s = 2mT, — 1.80 GeV/c^ 



<"[2m^, 1.80 GeV/c^] = (523.5 + 1.1 (stat.) + 1.0 (sys.) + 2.6 (Br.) + 2.5 (isospin)) x 10 



-10 



where the first error is statistical, the second is the experimental systematic (Table IXll) . the 
third comes from the uncertainties on the branching fractions (Table |X|) . and the fourth is 
from isospin-violating corrections (Table |XTT|) . 
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